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Abstract 

We introduce a new class of two-dimensional cellular automata 
with a bootstrap percolation-like dynamics. Each site can be either 
empty or occupied by a single particle and the dynamics follows a 
deterministic updating rule at discrete times which allows only emp- 
tying sites. We prove that the threshold density pc for convergence to 
a completely empty configuration is non trivial, < pc < contrary 
to standard bootstrap percolation. Furthermore we prove that in the 
subcritical regime, p < pc, emptying always occurs exponentially fast 
and that pc coincides with the critical density for two-dimensional ori- 
ented site percolation on T?. This is known to occur also for some 
cellular automata with oriented rules for which the transition is con- 
tinuous in the value of the asymptotic density and the crossover length 
determining finite size effects diverges as a power law when the criti- 
cal density is approached from below. Instead for our model we prove 
that the transition is discontinuous and at the same time the crossover 
length diverges faster than any power law. The proofs of the discon- 
tinuity and the lower bound on the crossover length use a conjecture 
on the critical behaviour for oriented percolation. The latter is sup- 
ported by several numerical simulations and by analytical (though non 
rigorous) works through renormalization techniques. Finally, we will 
discuss why, due to the peculiar mixed critical/first order character of 
this transition, the model is particularly relevant to study glassy and 



*Laboratoire de Probabilitcs et Modeles Aleatoires CNRS UMR 7599 Univ. Paris VI- 
VII 4,Pl.Jussieu F-75252 Paris Ccdcx 05, FRANCE; c-mail ctoninel@ccr.jussieu.fr 

^Service de Physique Theoriquc, CEA/Saclay-Orme des Merisiers, F-91191 Gif-sur- 
Yvette Cedex, FRANCE; email: giulio.biroli@cea.fr 



1 



jamming transitions. Indeed, we will show that it leads to a dynamical 
glass transition for a Kinetically Constrained Spin Model. Most of the 
results that we present are the rigorous proofs of physical arguments 
developed in a joint work with D.S. Fisher. 
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1 Introduction 



We introduce a new class of two-dimensional cellular automata, i.e. sys- 
tems of particles on 7? with the constraint that on each site there is at most 
one particle at a given time. A configuration at time t is therefore defined by 
giving for each x G the occupation variable ?7t(x) G {0, 1} representing an 
empty or occupied site, respectively. At time t = sites are independently 
occupied with probability p and empty with probability 1 — p. Dynamics is 
given by a deterministic updating rule at discrete times with the following 
properties: it allows only emptying sites; it is local in time and space, namely 
?7f+i is completely determined by r]t and ?7t+i(x) depends only on the value 
of r\t on a finite set of sites around x. 

We will be primarily interested in the configuration which is reached in 
the infinite time limit. We will first identify the value of the critical density 
Pc, namely the supremum over the initial densities which lead almost surely to 
an empty configuration. In particular we will prove that pc = where 
is the critical probability for oriented site percolation on 1? . Furthermore, we 
will analyze the speed at which the system is emptied in the subcritical regime 
and prove that emptying always occurs exponentially fast for p < 'p'^^ . Then, 
we will determine upper and lower bounds for the crossover length below 
which finite size effects are relevant when p Z' p^. These bounds establish 
that the crossover length diverges as the critical density is approached from 
below and divergence is faster than power law. Finally, we will analyze the 
behaviour around criticality of the final density of occupied sites, poo. We 
will prove that the transition is discontinuous: Poo(p) is zero if p < pc and 
Poo(Pc) > 0. We underline that both discontinuity and the lower bound on 
the crossover length are proved modulo a conjecture (Conjecture 13.11) for 
the critical behaviour of oriented site percolation (actually, for the proof of 
discontinuity we will only need a milder version of Conjecture 13.11 which is 
stated as Conjecture 13. 2p . This conjecture states a property which is due to 
the anisotropic character of oriented percolation and it is widely accepted 
in physical literature, where it has been verified both by analytical works 
through renormalization techniques and numerical simulations. However, we 
are not aware of a rigorous mathematical proof. 

One of the main interests of this new model relies on the peculiar feature 
of its transition: there is a diverging lengthscale as for standard continuous 
critical transitions, but at the same time the density of the final cluster 
Poo{p)i that plays the role of the order parameter, is discontinuous. This 
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discontinuous/critical character, to our knowledge, has never been found so 
far in any cellular automaton or in other type of phase transitions for short 
range finite dimensional lattices 0. 

Among the most studied cellular automata we recall bootstrap percola- 
tion [1] and oriented cellular automata [15]. In bootstrap percolation the 
updating rule is defined as follow^: a site can be emptied only if the num- 
ber of its occupied nearest neighbours is smaller than a threshold, /. In this 
case, it has been proved [1, 15] that on Z'^, the critical density is either 1 or 
depending on /: Pc = for I < d, pc = I for I > d. On the other hand, 
oriented cellular automata on Z'^ are defined as follows: site x can be emptied 
only if (x + ei, . . . , X + Crf) are all empty, where are the coordinate unit 
vectors. In this case it has been proven [15] that the critical density coincides 
with the critical probability for oriented site percolation and the transition 
is continuous, namely Poo(Pc) = 0. 

Our model shows a behaviour that is different from both bootstrap and ori- 
ented cellular automata, since the transition occurs at a finite density and it 
is discontinuous. Models with such a critical/first order transition have long 
been quested in physical literature since they are considered to be relevant 
for the study of the liquid/glass and more general jamming transitions. In 
the last section we will discuss the behavior of a Kinetically Constrained Spin 
Model, the so called Spiral Model (SM) [18], which has a stochastic evolu- 
tion with dynamical rules related to those of our cellular automaton. We 
will show that the present results for the cellular automaton imply that SM 
has a dynamical transition with the basic properties expected for glass and 
jamming transitions. For a more detailed discussion of the physical problem 
we refer to our joint work with D.S.Fisher [17, 18]. Most of the results that 
we present are the rigorous proofs of physical arguments developed in [17,18] 
for several jamming percolation models. The cellular automaton we consider 
in this paper (and the related Spiral Model [18]) is one of the simplest in 
this class. Originally, in [17], we focused on the so called Knight models for 
which some of our physical arguments cannot be turned into rigorous ones 

^On the other hand such type of transition is found in some problems on non- finite 
dimensional lattices, e.g. the k-core problem [13] or bootstrap percolation on random 
graphs [6] . It has also been established for long-range systems in one dimension [2] . 

^Note that the model is usually defined in this way in physical literature, while in 
mathematical literature the role of empty and filled sites is exchanged: dynamics allows 
only filling sites and a site can be filled only if the number of its neighbours is greater than 
I. The same is true for oriented models defined below. 
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as pointed out in [19]. In particular our original conjecture [17] that the 
transition for Knights occurs at p"^^ should not be correct [18,19]. However, 
as discussed in [18], numerical simulations suggest that the physical behavior 
of the Knight models around its transition (which is located before p^^) is 
the same of SM. 



2 Setting and notation 



2.1 The model 

The model is defined on the 2-dimensional square lattice, 1?. We denote 
by ci and 62 the coordinate unit vectors, by x, z the sites of 1? and by 
\x — y\ the Euclidean distance between x and y. The configuration space is 
Vt = {0,1}^^, i.e. any configuration 77 G f2 is a collection {r]{x)}x£z^2, with 
ri{x) G (0, 1), where and 1 represent an empty or occupied site, respectively. 
At time t = the system is started from a configuration 770 G chosen at 
random according to Bernoulli product measure fi^, namely ?7o(x) are i.i.d. 
variables and fi''{'r]o{x) = 1) = p. Therefore p will be called the initial density. 
The evolution is then given by a deterministic process at discrete time steps 
t = 0,1,2,... and the configuration at time t, rjt, is completely determined 
by the configuration at time t — 1 according to the updating rule 

Vt = Tvt-i (2.1) 
with the evolution operator T : Q ^ Q defined as 

if ri{x) = 

Tr]{x) ■■= { if rilx) = 1 and r/ G A (2.2) 
if ri{x) = 1 and r] ^ Ax 



with 

where 



A := (£,^^-^^ n ^-^^) (2.3) 



cNE-SW _ V in; 
cNW-SE _ 1; 111; 
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and, for any A C Z^, we denote by Va the event that all sites in A are empty, 
Va ■ iv ■ V{x) = \/x e A), and the sets NE^, SW^, NW^ and SE^ are 
defined as 

NE^ := {x + 62, X + ei + 62), 

SWx := (x - 62, X - ei - 62), 

A^H4 := (x - ei, X - ei + 62), 

^i^a; := (x + ei, X + ei - 62). 

In words the updating rule defined by (12.11) and (12.21) can be described as 
follows. Let the North-East (NEJ, South-West (SWJ, North-West (NW^) 
and South-East (SE^) neighbours of x be the couples depicted in Figure 
dJ If X is empty at time t — 1, it will be also empty at time t (and at any 
subsequent time). Otherwise, if x is occupied at time t — 1, it will be empty 
at time t if and only if at time t — 1 the following local constraint is satisfied: 
both its North-East or both its South-West neighbours are empty and both 
its North-West or both its South-East neighbours are empty too. See Figure 
(U^a) (Figure [T](b) ) for an example in which the constraint is (is not) satisfied. 
As it will become clear in the proofs of Theorems 13.41 and 13.51 the fact that 
in order to empty x we necessarily have to satisfy a requirement in the NE- 
SW and an (independent) requirement in the NW-SE direction is the key 
ingredient which makes the behaviour of this model quantitatively different 
from the oriented cellular automata in [15]. 

One can also give the following alternative equivalent definition of the 
dynamics. Let 2^ be the collection of the four subsets of each containing 
two adjacent couples of the above defined neighbours of x, namely 

:= {NE^ U SE^; SE^ U SW^; SW-^ U NW^; NW^ U NE^}. 

With this notation it is immediate to verify that definition (12. 3p is equivalent 
to requiring that at least one of the sets A G X,; is completely empty, namely 

The following properties can be readily verified. The dynamics is attrac- 
tive with respect to the partial order 77^ -< rf if ?7^(x) < rf{x) V x G Z^. 
Attractiveness here means that if we start the process from two different 
configurations rjl and ?7o with rjl < ijq at each subsequent time the partial 
order will be preserved. The updating rule is short-range, indeed Ax depends 
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only on the value of rj on the first and second neighbours of x. Furthermore 
it is both invariant under translations and under rotations of 90 degrees 
(and multiples). Indeed, if for all y G we define the translation operator 
Ty : fl —>■ Q as {Tyri)z = 'rjy+z, it is immediate to verify that Ty-r] G {^)Ax+y if 
and only if r/ G {^)Ax- On the other hand, if we let /„9o :7} 1? be the 
operator which acts as j-oii^{x\e\ + X2t2) := —X2t\ + X\e2 and we define the 
rotation operator T^go : f2 ^ ^2 which acts as (7?.9o?7)(x) = ?7(/_9o(x)), it is 
immediate to verify that 'R-^^r] G (^)^/_go(a;) if and only if 77 G {^Ax- 

For future purposes it is also useful to define the model on a finite volume 
A C Z^, i.e. to define an evolution operator Ta : ^a where f^A is the 

configuration space fi^ := {0, l}'^. A natural way to do this is to fix a 
configuration uj G and to consider the evolution operator Ta^o; with 

fixed boundary condition tu, i.e. for each x G A we let 

Tf,^^ri{x):=T{ri-uj){x) (2.4) 

where 77 ■ G is the configuration which equals 77 inside A and uj outside. 
Note that Ta,^^ depends only on the value of lo on the sites y ^'L\K such 
that y G [NEx U SEx U SEx U NWx) for at least one x G A. Note also that 
the configuration reached under T^^.a after |A| steps is stationary, namely 

^> = (2.5) 

for any n > (this trivially follows from the fact that we are evolving deter- 
ministically on a finite region and that only emptying of sites is allowed). A 
choice which we will often consider is the case of filled boundary conditions, 
namely u;(x) = 1 for all x G Z\A, and we will denote by T/ the corresponding 
evolution operator. 



2.2 Main issues 

Before presenting our results, let us informally introduce the main issues 
that we will address. We underline once more that these are akin to those 
examined in previous works for bootstrap percolation and oriented models 
[5,15]. However, the answers will be qualitatively different. 

• We will determine the critical density pc such that, a.s. with respect 
to the initial distribution /i'', if p < all the lattice gets eventually 
emptied under the updating rule, while for p > Pc this does not occur. 
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Figure 1: Site x and the four couples of its North-East (NE), South-East 
(SE), North- West (NW) and South- West (SW) neighbours, (a) If x is occu- 
pied it will be empty at next time step. Indeed its NE and NW neighbours 
are all empty, thus t] G Vne C S'^^-^^ and rj e Vnw C. g^'^-^^ , Therefore 
Tj G Ax- (b) If X is occupied it will remains occupied at next time step. 
Indeed neither the NE nor the SW neighbours are completely empty, thus 
7] ^ g^E-sw ^T^^ therefore 77 ^ Ax- 



The precise definition of pc follows. Consider on {0, 1} the discrete 
topology and on VL the Borel a-algebra S. Let M. be the set of measures 
on (f2, S) and /i^ be the evoluted of the initial distribution /Xq = 
according to the above deterministic rules. Due to attractiveness it 
is immediate to conclude that /if converges weakly to a probability 
distribution G M.- Following [15] we can indeed define a partial 
order among /i, z/ G 7W as /i < z/ if / fijf) dnirf) < J f{ri)di'{ri), V/ : 
r2 — s> R and / increasing. The fact that O's are stable implies that /Iq > 
y^i > • • • . This, together with the compactness of Q and Ai, assures 
the weak convergence of /xf to a probabihty distribution /i^ G ^A . We 
can therefore define the critical density pc as 

Pc := sup (p : Poo (p) = 0) (2.6) 

where poo(p), henceforth referred to as the final density, is defined as 

Poo(p):=/iSo(^(0)) (2.7) 



• We will analyze the speed at which the system is emptied in the sub- 
critical regime, p < pc- Let tshe the first time at which the origin gets 
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emptied 

tE ■■= inf (t > : ritiO) = 0). (2.8) 

Following notation in [15] we let -Pp(-) be the probability measure on 
({0, i}^^{o,i,2,...}^ „here is the Borel a-algebra on {0, i}^'.{o,i,2,...}. 
With this notation we define the speed 7 as 

7(p) := sup (/5 > : 3C < 00 s.t. Pp{tE > t) < Ce"^*) (2.9) 

and the corresponding critical density as 

= sup (p : 7(p) > 0) (2.10) 

It is immediate from above definitions to check that pc < pc- We will 
prove that the equality is verified, namely emptying always occurs ex- 
ponentially fast in the subcritical regime p < pc- 

• We will analyze the final density at criticality and establish that the 
transition is discontinuous (poo(Pc) > 0). 

• We will analyze the finite size scaling. Let K2L C 7? and C Al 
be two square regions centered around the origin and of linear size 2L 
and iv/2, respectively. We denote by ri'^{ri) the stationary configura- 
tion which is reached after (2L)^ steps when we evolve from rj with 
filled boundary conditions on A2L, ^^{rf) ■= {tJ^^^Y^^ r]. Finally, we let 
E{L,p) be the probability that Al/2 is empty in r]'^ 

E{L, p) = fi'iv'ix) = Vx G Al/2) (2.11) 

As we shall show lim^^^oo E{L, p) = 1 for p < pc and limpyp^ E{L, p) 7^ 
1 when L is kept fixed to any finite value. We will therefore study the 
scaling as p pc of the crossover length H(p) defined as 

S(p):=inf(L:E(L,p)>l/2) (2.12) 

Another possible definition of the crossover length would have corre- 
sponded to defining -E(L, p) as the probability that the origin is empty 
in ?7*. This requirement, less stringent than the previous one, leads to 
the same results for S(p) at leading order, see Section 6. 
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3 Results 



Let us first recall some definitions and results for oriented site percolation 
on 1? which we will be used in the following. We say that (xi,X2, • • -Xn) is 
an oriented path in if Xj+i — Xj G (61,62) and, given a configuration 
r] G {0, 1}^^, we let x — > ?/ if there exists an oriented path connecting x and 
y (i.e. with xi = x and Xn = y) such that all its sites are occupied (rjixi) = 1 
for all z = 1, ... ,77.) . For a given x G Z^, we define its occupied oriented 
cluster to be the random set 

C^''iv)--=iyeZ':x^y) (3.13) 

Finally, for each n, we define the random set as 

F" := (y G : X ^ y and 3m s.t. y — x = mci + 77,62) (3-14) 

With the above notation the percolation probability a{p)^^ is defined as 
a{p)^^ := /i^(?7 : | C^^ \ = 00). As it has been proven, see [7], a{p)'-^^ is zero 
at small p and strictly positive at high enough p: the system undergoes a 
phase transition. The critical density, defined as p'^^ := inf(p : a(p)'^^ > 0), 
has been proven to be non trivial, < p'^^ < 1, see [7] for some upper 
and lower bounds. We also recall that extensive numerical simulations lead 
to p^^ ~ 0.705489(4) [11]. Furthermore the transition is continuous in the 
percolation probability, namely a{p^^) = [8] and in the subcritical regime 
an exponential bound has been proven [7]: at any p < p^^ there exists 
^op{p) < 00 such that for n ^ 00 the following holds 

/i^(F^ ^ 0) < e""/«o^. (3.15) 

Finally, we recall the conjecture for the critical behavior. Let Aa^b be a 
rectangular region with two sides (i9-Ri and 9i?2) of length a parallel to 
61 + 62, and two sides {dR^ and dRi) of length b parallel to —61 + 62. Let also 
f^Nz Bernoulli measure on A7v,Ar^ conditioned by having both sides 

of length A^, dTZi and dTZ2, completely empty. The following properties are 
expected to hold for the probability of finding an occupied oriented path 
crossing the rectangle in the direction parallel to 61 + 62, i.e. connecting the 
two non empty borders dR^ and 9i?4 
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Conjecture 3.1. There exists z, c^p, Cqp and ^(p) with < 2 < 1, < 
Cqp < 1, < Cqp < 1, ^(p) < 00 for p < p^^ and liiiip^pop ^ = 00 s.t. 

O P 

lim ( 3 a; G 97^3 and y G dTZ^ s.t. x ^ y) = Cqp (3.16) 

L— >oo ' 

lim /if X G dTZ-i and y G dTZ^ s.t. x ^ y) = Cqp (3-17) 

The above conjecture is given for granted in physical literature, where 
it has been verified both by numerical simulations and by analytical works 
based on renormalization techniques (see [11] for a review). The physical 
arguments supporting this conjecture are based on finite size scaling and the 
anisotropy of oriented percolation that gives rise to two different correlation 
lengths in the parallel and perpendicular direction w.r.t. the orientation of 
the lattice, namely in the ei + 62 and — ei + 62 directions. This explains why 
in finite size effects an anisotropy critical exponent z emerges such that at 
P = P*^^ the probability of finding a spanning cluster on a system of finite 
size Lx converges to a constant which is bounded away from zero and one 
when L 00. To our knowledge, a rigorous proof of Conjecture 13. II has not 
yet been provided. However, in [7] it has been proven that the opening edge 
of the percolating cluster is zero at criticality, which implies anisotropy of the 
percolating clusters. Our results, unless where explicitly stated (Theorem 13.41 
and 13.51 (ii)), do not rely on the above Conjecture. 

Finally, we recall that the minimal C,op and ^ for which fl3.15p and (13.161) 
hold, namely the parallel correlation length, is expected to diverges when 
P y P?^ a,s (p — p^^)" with a ^ 1.73. The lowest value of z for which 
the result of Conjecture 13.11 is expected to hold, i.e. the anisotropy critical 
exponent, is z ~ 0.63. 

Before stating our results, let us give a milder version of Conjecture 13.11 
which will be sufficient to prove discontinuity of the transition (the stronger 
version 13.11 will be used only to prove the upper bound for the correlation 
length). 

Fix io > and consider the sequence of increasing rectangles TZi : = 
^ii,i/i2£i with £i = 2£j_i and denote the two short sides parallel to the —61+62 
direction by dR^ and dR^. Let Si be the event that TZi is crossed in the 
parallel direction, namely 

Si:={r]: 3x E dUl and y G dUl s.t. x ^ y) (3.18) 
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Conjecture 3.2. X]j=i I log(/i''(^,))l <^ at p = 



The fact that 13.11 imphes 13.21 follows immediately by cutting TZi into 
0{il~^) slices of size ii x £f and using fl3.16p to bound the probability that 
each slice is not spanned by a cluster. 

We are now ready to state our results. We have proved that the critical 
densities defined in (12. 6 p and fl2.10p are equal and furthermore they coincide 
with the critical probability for oriented site percolation (and therefore also 
of oriented cellular automata [15]) on Z^, namely 

Theorem 3.3. pc = Pc = P^^ 

However the critical properties are different from oriented percolation: 
the transition is here discontinuous in the final density and the crossover 
length diverges faster than any power law at criticality. More precisely 

Theorem 3.4. // Conjecture ] 3.^ holds. Poo(Pc) > 0. 

Theorem 3.5. i) limpyp^ ^op(p)~^^'' logS(p) = for any e > 0. 

a) If Coniecture \3.1\ holds, S(p) > ci^ exp[c2 ^(p)^~^] 
with ci = 1/(2^2) and ca = | log(l - c^op)\/'^ 

where ^op{p) is the smallest constant for which (13.150 holds and ^ are 
the smallest constant which satisfy (I3.17p . Note that, if the conjectured 
power law behavior for ^ and ^qp holds then our bounds imply a faster than 
power law divergence for H. This property, as well as the discontinuity of 
the final density, makes the character of this transition completely different 
from the one of oriented percolation and oriented cellular automata. 



4 Critical density: proof of Theorem 13.3 



Proof of Theorem \3.3\ The proof follows from the inequality p^ < Pc (which 
can be readily verified from definition 12.61 and I2.10p and the following Lemma 
OandlS □ 

Lemma 4.1. pc < P^^ 

Lemma 4.2. 7(p) > for p < p^^ ■ Therefore pc > p^^ ■ 
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4.1 Upper bound for pc,: proof of Lemma 14JJ 



In order to establish an upper bound for pc we first identify a set of config- 
urations in which the origin is occupied and it can be never emptied at any 
finite time because it belongs to a proper infinite cluster of occupied sites. In 
this case we will say that the origin is frozen. Then we prove that a cluster 
which makes the origin frozen exists with finite probability under the initial 
distribution /Xq = for p > . This follows from the fact that the origin 
can be frozen via two infinite independent clusters which, under a proper 
geometrical transformation, can be put into a one to one correspondence 
with infinite occupied clusters of oriented percolation. We stress that these 
clusters, which are sufficient to prove the desired upper bound for the critical 
density, are not the only possible clusters which can freeze the origin, as will 
become clear in the proof of Theorem 13.41 

Let us start by introducing some additional notation. 
We say that (a;i, X2i ■ ■ ■ Xn) is a North-East (NE) path in if Xj+i G NE,;^.. 

for alH G (1, . . . , n — 1), and we let x — > y if there is a NE path of sites 
connecting x and y (i.e. with xi = x and x„ = y) such that each site in 
the path is occupied {rjlxi) = 1 for all z = 1, . . . ,n). Also, we define the 
North-East occupied cluster of site x to be the random set 

C!:^:={yeZ^:x^y) 
and, for each n, we also define the random set F^"^'" as 



r^^'" := (y G : X ^ y and 3m s.t. y - x = mci + ne2) (4.19) 

We make analogous definitions for the South-West , North-West and 
South-East paths, the correspondent occupied clusters , and C^^ 

and for the random sets F;^^'", T;^^'" and T;^^'" . Note that if x is empty its 
occupied cluster in all the directions is empty, instead if x is occupied each 
occupied cluster contains at least x. In Figure [2] we depict as an example the 
North-East (inside the dashed line) and South-West (inside the continuous 
line) occupied clusters of a given occupied site x. With this notation we 
define 

^NE-SW . I ^NE I ^ ^ I ^SW I ^ ^) (4 20) 

and it is immediate to verify that the origin is frozen on any configuration 
T] G , namely 
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Figure 2: Sites inside the continuous (dashed) hne form the South- West 
(North-East) occupied cluster for the origin, 0. It is immediate to check that 
in order to empty we have (at least) to destroy either the NE occupied 
path connecting to |/ or the SW occupied path connecting to w (sites 
indicated by the dotted line). This requires a number of steps which is at 
least the minimum of the lengths of these two paths (each path have to be 
emptied sequentially from its external border unless the other one has been 
already emptied). 



Lemma 4.3. p^{p) > fiP{J^^^-^^) 

Proof. The result follows immediately once we prove that, for any given 
r > 0, the inequality Ie > t holds, where Ie is the first time at which the 
origin gets emptied, see definition (12. 8p . Since for hypothesis \Cq^\ = oo and 
\Cq^\ = oo, there exist y and w such that y G C^^, w G Cq^, the length 
Ly of the minimal NE occupied path xi = 0, . . . ,XLy = y connecting to 
y verifies Ly > t + 1 and the length Ly^ of the minimal SW occupied path 
xi = 0, . . . ,xl^ = w connecting to w verifies > r + 1. Let X2 be the 
site which is emptied first among X2 and X2. It is immediate to verify (see 
Figure [2]) that > inf(t > : ?7t(x2) = 0) > min(L^„,Lj^) — 1 > r, which 
concludes the proof. □ 

By definition of NE and SW neighbours, it is easy to verify that (except 
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for X itself) there do not exist sites that can be connected to x both by a NE 
and a SW path, thus it follows immediately that 

^P^:pNE^sw^ > /i-d I = oo) /i^'d Cl^ I = oo) (4.21) 

We will now prove that the probabilities of such infinite NE or SW occupied 
clusters can be rewritten in terms of the probability of infinite clusters for 
oriented percolation, namely 

Lemma 4.4. C^^ | = cx)) = C^^ | = oo) = yu(|Co°^| = oo) 

Proof. Let Vi = 61+62 and f 2 = 62, each x G can be written in a unique 
way as X = mvi + nv2. We can therefore define the operator R'^^ : f2 — *• O 
which acts as {R^^ri){mei + 71,62) = ri{mvi + nv2)- It is immediate to verify 
that nP{r]) = fiP{R^'^7]) and that |C^(r/)| = \C^^{R^^r])\ for any r/. Thus 
^p(\ C^E I _ ^) _ ii{\C^P{r])\ = 00). The result C^^ | = 00) = 
/i(|C^^| = 00) can be proved analogously. □ 

We are now ready to conclude the proof of the upper bound for p^. 
Proof of Lemma \4.1\ The result follows from Lemma 14.31 equation fl4.2ip . 
Lemma |4]4] and the definition of p^^ which implies /i(|C^^(?7)| = 00) > for 
p>p?r □ 



4.2 Lower bound for pc'. proof of Lemma 14.21 

The central result of this section is Lemma 14171 This contains a lower bound 
for the probability that a certain finite region can be emptied (except for some 
special sets at its corners) when evolution occurs with fixed filled boundary 
conditions and p < . Since this lower bound can be made arbitrarily near 
to one provided the size of the region is taken sufficiently large, the result 
Pc > P^^ will easily follow (Corollary 14. 8p . Some additional work involving 
a renormalization technique in the same spirit of the one used for bootstrap 
percolation in [15] will be used to prove the stronger result > p^^ (Lemma 

USD. 

Let Sa be a segment of length a with left vertex in the origin, 

Sa := U^.-i^6i 

and Tlafi be the quadrangular region with two sides parallel to the 61 direction 
and two sides parallel to the 61 + 62 direction which is obtained by shifting 
h times Sa of 61 + 62 depicted in Figure [3|, namely 

T^a,b ■■= Uto [Sa + ^(ei + 62)] 
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Figure 3: (a) The region 7la,b (here a = 7, b = 6) and (inside the continuous 
hue) the boundary region d7la,b which guarantees that we can empty 7la,b 
(Proposition 14.51) . (b) The alternative choice of boundary conditions which 
is described in Remark 14.61 

where for each x G and A G 1^"^ we let x + A G Z"^ he x + A := {y : y = 
x + z with z & A). As it is immediate to verify, if we impose empty boundary 
conditions on the first external segment parallel to the bottom border and 
on the first two segments parallel to the right border (empty sites inside the 
continuous line in[3](a)), 7la,b is completely emptied in (at most) |7^a,6| = o,b 
steps. More precisely, if we define the bottom right border as 

dna,b ■■= {Sa-e2-ei)U[aei + {b - l)(ei + e2)]uti[^2 + ae^ + {i - 2)(ei + es)] 

and we recall that, for each A C Z^, Va is the set of configurations which are 
empty on all sites in A, the following holds 

Proposition 4.5. If rj e Vgn^ t, then T"''?] G V-r,,^. 

Proof. Starting from the bottom right corner of TZa,b we can erase all 
particles in Sa from right to left, thanks to the fact that their SE and SW 
neighbours are empty. Then we can erase all particles in Sa + e2 + ei starting 
again from the rightmost site and so on until emptying the whole region. □ 

Remark 4.6. Analogously, it is easy to verify that an alternative choice of 
boundary conditions which guarantees that we can empty TZa,b is to impose 
empty sites on the external segment parallel to the top border and on to the 
two external segments parallel to the left border (see Figure\^b)), namely on 

{Sa + L{ei + 62)) U (-ei) Uti [^2 - 2ei + z(ei + 62)] 

With a slight abuse of notation we denote the region with the shape of 
TZl,l which is centered around the origin and its corresponding border by 
TZl and dRi. 

7^L := 7^L,L - L/2e^ - L/2e2 
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ORl := dnL,L - L/2e, - L/2e2 

namely TZl is the region delimited by vertexes A,B,C,D in Figure [7] (here 
and in the following, without lack of generality, we choose L such that L/4 
is integer). Then we denote by TZl the region inside the bold dashed line in 
Figure U\ which is obtained from TZ^ by subtracting at the bottom left and 
top right corners two regions, TZu and TZtr, which have the shape of TZl/a, 
namely 

■JZl :=7^L\(7^MU7^^,) 

with 

TZtr := 'R-L/i,L/i + L/A{ei + 62) 
Tibi ■= TIl/4,l/4 - Lei - L/2e2 

Let ?7 be a configuration on TZl and denote by 77^ the stationary configuration 
reached upon evolving 77 with fixed filled boundary conditions on 7^^, 77* : = 
(T^^) 77. We say that r] is good if TZl is completely empty on r]'^, and we 
denote by the set of good configurations, i.e. 

:= iv e : v' e V^J. 

Given a configuration r] on we denote by rjn^ its restriction to TZl. If 7/7^^ 
is good, then the evolution on the infinite lattice also empties in at most 
steps the region TZl, namely r] E V^^, as can be easily proved by using 
attractiveness of the dynamics. 

The following holds on the probability that a region is good 

Lemma 4.7. For any p < p^^ and for any e > there exists L{p, e) < 00 
such that for L = L{p, e) 

fiP{G^) > 1 - e 

We postpone the proof of this main Lemma 14.71 and derive its conse- 
quences for pc and Pc- 

Corollary 4.8. pc > p^^ , thus pc = p'^^ ■ 

Proof. Since the origin belongs to TZl, ii TZl is good then the origin is 
certainly empty at time L^, which implies 

/i''(T^'r/(0) = 0) > p'^iGL) (4.22) 
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This, together with the definition 02.71) and the result of Lemma 14. 7[ guar- 
antees that if p < for any given e > we can choose L > L{p, e) such 
that 

1 - f^Livm > /""(T^'^lO) = 0) > 1 - e (4.23) 

Thus for any e we have < /i^(?7(0)) < e and therefore /i^^(?7(0)) = 0, which 
imphes pc > p^^ (recall definition fl2.6p for pc). The identification of pc with 
immediately follows from the latter result and Lemma 14.11 □ 

In order to prove the stronger result of Lemma [4.21 we now have to intro- 
duce a renormalization procedure in the same spirit of [15]. Fix an integer 
scale L and let 1?{L) = LI? . We consider a partition of 1? into disjoint 
regions 7^|^ := 7^^ + 2:, 2; G 7}{V). In the following we will refer to iT'iV) as 
the renormalized lattice and, given configuration rj G ^2^2 , we will say that a 
site z G Z^(L) is good if the configuration 777^^ restricted to the tile 7?.|^ C 
corresponding to z is good. Note that the events that two different sites z 
and z' are good are independent. 

Let z be site of the renormalized lattice. If its South, SouthEast and East 
neighbours, i.e. z — 62, z + ei — 62 and z + ei, are good then after (at most) 
2|7^L| + |7^L/4| = ^233/16 steps the region corresponding to z on the original 
lattice, TZf^, is completely empty. More precisely 

Proposition 4.9. If rj^z+e^ G Gl, r/^^+ei-ej G Gl and G Gl then 

T^^'v e Vni for C = 33/16. 

Proof. From the definition of good configurations it follows that at time 
L2 all sites belonging to TZl^^' U ni"^^ U TZl^^'-^^ are empty, i.e. 

G Va with A := ni+^' U ni~'^ U ni+''-'^ (4.24) 

If we denote by TZtr^^ (J^li^^^) the top right (bottom left) region of linear size 
L/4 which belongs to TZl'^^ (7^^'') and by dTZt'^^ {dKt'') the correspond- 
ing bottom right borders, it is immediate to verify that fl4.24p implies that 
at time both d7l^~^^ and dTZl^'^^ are completely empty. This, together 
with Proposition 1131 implies that at time + (1/4)^ both TZ^''''' and TZl^''^ 
will be empty. The latter result, together with (14.241) . guarantees that the 
bottom right border of 7^|_, 97?.|_, is empty at time + (L/4)^. By using 
again Prop osition 14 . 51 it follows that at time L^ + L^/16 + L^ the entire region 
TZf^ will also be empty. □ 
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Figure 4: The renormalized lattice: black and white circles stand for bad and 
good sites, respectively. Inside the continuous line we depict the S-SE-E bad 
cluster of the origin, C^. For this configuration the range of the bad cluster 
verifies = 7. 

We say that {xi,X2, ■ ■ ■ Xn) with Xi G 1?{L) is a South-SouthEast-East 
(S-SE-E) path if Xj+i G (ei, ei — 62, —62) for alH G (1, . . . , n — 1). For a given 

configuration 77 let x -—^ ^ y if there is a S-SE-E path of sites connecting 
X and y (i.e. with xi = x and Xn = y) such that each site in the path is not 
good. Also, we define the S-SE-E had cluster of the origin to be the random 
set 

C'':= (yeZ\L):0^-^''y) 
Finally, we define the range of the bad cluster of the origin, i?'', as 

:= sup(A; : n {Dk U Dk+i) ^ 0) 

where 

Dk = {y ■ y = a{ei + 62) - ke2 with a G (0, . . . , /c)) 

(and we let sup(0) = —00). See figure H] for an example of C'^ and R^. 
Proof of Lemma \4.2[ 

Let r be the renormalized time defined by the following relation 

t(r) := 33/16LV (4.25) 

For a given configuration 77, we denote by and R^ the cluster and range 
of the origin at the renormalized time r, i.e. for the configuration ?7t(T-)- Let 
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i?Q = k. At time r = 1 Proposition 14.91 guarantees that all the tiles of the 
original lattice corresponding to renormalized sites in and which belong 
to the cluster of the origin are empty. Indeed each of these sites has its South, 
SouthEast and East neighbours which are good at time zero since they are 
in Dk (if they were bad they would also belong to the cluster of the origin 
and this would imply > k + 1 in contrast with our assumption Rq = k). 
This implies therefore R\ < k — 1 = Rq — 1. The same procedure can be 
applied at each subsequent (renormalized) time step, yielding 

< - 1 (4.26) 

and finally 

Rl < 0. (4.27) 

The latter equation implies that the South, the SouthEast and the East 
neighbours of the origin are good when the renormalized time coincides with 
the range of the bad cluster of the origin at time zero. Therefore, using 
again Proposition 14.91 we get that the origin is certainly empty at time 
t = 33/16L^(i?Q + 1) and therefore the first time at which the origin gets 
emptied, Ie, verifies 

Pp{tE >t) < ^i''{R^ > tl6/33L~^ - 1) (4.28) 

We can now use a Peierls type estimate to evaluate the probability that the 
range of the origin is larger than a certain value 

oo 

fi^iR' > s - 1) < ^(1 - fiP{GL)y^^3' (4.29) 

l=s 

If p < p^^ and we choose e < 1/3, Lemma (14. 7p guarantees the existence of 
an L(p, e) such that if we let L = L{p, e) we have a := {1—p''{Gl))3 < 3e < 1. 
Therefore we obtain and exponential decrease for the above probability, 

pf{R^ > s - 1) < (1 -a)^iexp(-s|log(a)|). (4.30) 

This, together with (14.281) . allows to conclude 

Pp{tE>t)<Cexp{-t(3{p)) (4.31) 

with j3{p) = I log(Q!)|16/33L(p, 1/3)"^. Thus the the speed 7(p) at which 
the lattice is emptied (see definition fl2.9p ) satisfies 7(p) > /3(p) > at any 
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p < jf^^ and we conclude that > p^-^. Note that, as a byproduct, we have 
also derived a lower bound on the velocity in terms of the crossover length 
of oriented percolation (via expression fl4.4ip for L{p,e)). □ 

We are now left with the proof of the main Lemma 14.71 which will be 
achieved in several steps. Let us start by proving some results on the sufficient 
conditions which allow to enlarge proper empty regions. 

Let Q^^-^^ and Q^^-^^ be the two quadrangular regions inside the 
continuous lines of Figure [5l 

qNW-ne ^ ui,(52L-2(^-i) -{L-Z + l)ei + (z - 1)62) 
qSW-se ^ ui,(52L-2(.-i) -{L-Z + l)ei - ze^) 
and Ol be the octagon centered around the origin formed by their union 

Ol := Q^^-^^ u Ql^-^^. 

If Ol is empty a sufficient condition in order to expand the empty region of 
one step, i.e. to empty the region is that the four key sites 

Kr ■■=-iL+l)eu K^^:=e, + Le2 

Ki"^ := -{L + l)e„ Kf^ := {L + l)e, - 62 



are all empty, namely 

Proposition 4.10. // G Vo^ and 7]{K^^) = r]{K^^) = r]{Kl^ 



ri{Kl^) = 0, then T^t] E Vo, 



Proof. From Proposition 14.51 it is immediate to verify that we can empty 
all sites of the form + ^ + ^2) with 1 < z < L + 1. The same 

occurs for all sites Kl^ ~ ~ + G2) with again 1 < z < L + 1, as can 
be immediately be verified by using Remark 14.61 Then it is easy to verify 
that we can subsequently empty all sites of the form K^^ + ~ " ^2) 
[Kl^ + (z - l)(-ei + 62)) with 1 < z < L + 1 since they all have the NW 
and SW {NE and SE) neighbours which are empty. □ 

We will now prove that, if Ol is empty, a sufficient condition in order to 
guarantee that is empty after steps is that its NE occupied cluster 



does not survive after L steps, namely r^;^, is an empty set (and analogous 
results in the other directions). More precisely if we define the events 

:= iv ■■ r'^fw = 0), := iv ■■ 4- = 0) (4-32) 
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Figure 5: The octagon Ol composed by Q^^-^^ and Q^^z-Si? ^^^^ 
bottom regions inside the continuous hues, respectively) and the key external 
sites i^f^,i^f^,i^f^,i^f^. 

:= (r^ : Vf,,' = 0), JC^ := iv : ^^^7.' = 0), 
the following holds 

Lemma 4.11. (z) If r] e Vo^ and r] e /Cf^ then T^\{K^^) = 0. 
(zi) Ifr]e Vo^ and r] e /C^^ then T^'7]{K^^) = 0. 
(m) Ifr]e Vo^ and t] E Kf^ then T^'r]{Kl^) = 0. 
(iv) IfriE and r] G /Cf^ then T^\{Kl^) = 0. 

Proof. We will prove only result (i) , the proofs of the other results follow 
along the same lines. Let V be the square region of size L x L with vertexes 
(— ei, , + Le2, —ei + Le2) (region inside the dashed line in figure [5]). 
Recall the definition fl2.4|] and consider the evolution operator p restricted 
to V and with fixed boundary conditions uj{x) = for x G Ol and u;(x) = 1 
for X G \ Ol. For simplicity of notation we will call T such operator. It 
is immediate to verify that 

If r/ G Vo^ and ^^^{K^^) = 0, then T'^'r]{K^^) = 0. (4.33) 

We are therefore left with proving that the hypothesis in (i) imply those in 
fl4.33p . which we will do by contradiction. First we need to introduce some 
additional notation. 
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Let 77** be the stationary configuration reached under T after steps, 
T]" := T^^ (stationarity follows from fl2.5p and \V\ = L"^) and V* be the 
rectangular region V* : V U,^^q ne2 U^Iq {—nei — 62), we define the following 
random set 

B:={xeV: r]'{x) = 1; (4.34) 
ri^{y) = if y eV* and y = x + {n + m)ei + me2 with n > 1, m > 0; 
V^{y) = if ?/ G P* and y = x + {n + m)ei — (m + l)e2 with n > 1, m > 0) 

The following properties, whose proof is postponed, hold 

Proposition 4.12. If x E B, then 

(i) ri'^{x + 62) = 1 or rj^{x + ei + 62) = 1 (or both); 

Let b{x) be rightmost among these two sites which is occupied (in r)'^) 

(li) Ifb{x) E V, then b{x) E B. 

By using the above properties, we will conclude the proof of lemma 14.111 
by contradiction. Let us suppose that the left hand side of (14.331) does not 
hold, namely rj'^^K^^) = 1. Since 77 E Vol follows immediately that 
1<[NW g ^ Thus we can define a sequence {xi} with i E (1,...,L) with 
Xi = and Xi = b{xi-i) for i < L and it is immediate to verify that 

Xi E NEx-_^ and rj^ixi) = 1. Therefore, under the hypothesis ri^{K^^) = 1, 
we have identified a NE occupied path xi, . . . ,xl of L sites for Kf^^ such 
that xl — = mci + Le2, which implies F^^;^ ^ 0. Since this result 

contradicts the hypothesis of the Lemma we conclude that ri'^{K^^) 7^ 1, 
thus the condition in the left hand side of (14.331) is verified and the proof is 
concluded. □ 

We are therefore left with proving the properties of the random set of 
sites B. 

Proof of Proposition \4.1^ 

(i) follows immediately from the fact that the South-East neighbours of x, 
SEx = (x + ei, X + ei — 62), are both empty in rj^ {x E B implies x eV and 
therefore SE^ C V* and both these sites are empty thanks to the definition 
of B). Thus at least one of the two sites NE^ = {x + e2, x + 61 + 62) should be 
occupied, otherwise x would be empty at time s + 1 (which is forbidden by 
the stationarity of 77**). We let b{x) be the rightmost among these occupied 
sites. 

(ii) The first property defining B, rj^{b{x)) = 1, is satisfied by definition of 
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b{x). Let us prove that the second and third property are verified too. If 

b{x) = X + ei + 62, then 

b{x) + (n + m)ei + me2 = x + {n + m + l)ei + (m + l)e2, 

b{x) + (n + m)ei — (m + l)e2 = a; + (n + 1 + m)ei — me2 
and the properties defining B are immediately verified. On the other hand, 

if b{x) = X + 62, 

b{x) + (n + m)ei + me2 = x + ((n — 1) + (m + l))ei + (m + l)e2 

+ (n + m)ei — (m + l)e2 = x + ((n + 1) + (m — l))ei — ((m — 1) + l)e2. 

Thus if n > 2 (for all m > 0) the second property defining B for site b{x) 
follows from the fact that x E B and the same holds for the third property if 
m > 1 (for all n>l). The third property is also easily established if m = 0, 
since in this case b{x) + (n + m)ei — (m + l)e2 = b{x) + nei — 62 = x + nei, 
which is again empty since x & B. Some additional care is required to verify 
the only remaining case, i.e. the validity of second property in the case 
n = 1. Notice that since b{x) = x + 62, this implies r]^{x + 62 + 61) = 
(by definition b{x) is the rightmost occupied site in NE^). Thus the second 
property is verified in the case n = 1 and m = 0. Let us consider the case 
b{x) + {n + m)6i + m62 with n = 1 and m = 1. It is easily verified that this 
site is also empty in the stationary configuration 77* since its SE neighbours 
are of the kind x+ (l+m)6i+m62 or x + (2+m)6i+m62 (and therefore empty 
since X G B) , one of its SW neighbours is also of the form x+ (l + m)6i + m62 
and the other one is b{x) + (1 + (m — l))ei + (m — 1)62. The latter, for the 
choice m = 1, is b{x) + 61 which we verified to be empty. The same property 
can be verified by induction for all the other values m > and n = 1. □ 
From Proposition 14.101 and Lemma [4.111 we can now derive the following 
lower bound (which will be used to prove Lemma HTT]) on the probability that 
after L„ := 4(L + n)^ steps we have enlarged of (at least) n steps the empty 
octagonal region 

Lemma 4.13. If p < v'^^ 1 then 

n-l 
i=0 
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Proof. 



Case n = 1. 

From Proposition 14.101 and Lemma 14.111 it follows that 



where we used the fact that the events V©^, /C^^, IC^^, JCf^ and /Cf^ are 
independent and we used the equalities 

/i^l/Cr ) = /i^(/Cr) = f^'ilCD = /i^(/CD (4.35) 

which easily follow from symmetry properties. Along the same lines of the 
Proof of Lemma 14.41 it is now easy to prove that the probability of this event 
coincides with the analogous quantity for oriented percolation, namely 

li'iJC^'^) = ^'{Tf^t = 0) = y^^ir^w = 0). (4.36) 



Then the proof is completed by using the exponential bound fl3.15p . 
Case n = 2 

From Proposition 14.101 and Lemma 14.111 we get 



/i''(TSeVo,^,|r/GVoJ> (4.37) 



where the couples of events at size L and L + 1 are now not independent, 
H^ilC^^ n /C^li) ^ i2p{IC^'^)i2p{IC^^^) However, since all the events that we 
consider are of the form fl4.32l) and therefore non increasing with respect to 
the partial order r/ -< r/' if ?7(x) < ri'{x) Vx, we can apply FKG inequality [10] 
and get 



/i^(T% e Vo^^Jr/ G VoJ > (4.38) 

[/i^(/Cr )]' [f^'ilC^Z)]' > [1 - ^M'L/^op)t [1 - exp(-(L + l)/^op)t 

where again we used the symmetry properties fl4.35p . the mapping to oriented 
percolation (14.361) and the fact that p < p^^. 
Case n > 2. 

The proof follows the same lines of the case n = 2. □ 
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Figure 6: Inside the dashed hne we depict the region 7^ := {S2L — Lei)U{e2 + 
qNE-nw-j which is emptied provided Q^-^-^^ (region inside the continuous 
hne) is empty and the other hypothesis of Lemma 14.141 hold. The region 
inside the dashed-dotted hne is A2m,n, the rectangular region whose internal 
configuration is not involved in the hypothesis of Lemma 14.141 



In the proof of the main Lemma W77\ we will also need a condition which is 
sufficient to guarantee the expansion of an empty region of the type Q^^-^'^ 
to the larger region 7^ := {S2L — Lei) U (62 + q'^^-^'^) (region inside the 
continuous line of Figure [6]). In particular, we will need a sufficient condition 
which does not involve the configuration inside the rectangular region A2m,n 
with vertexes (L + 2)e2 — (m — l)ei, (L + 2)e2 + mei, — (m — l)ei + (L + 2 — 
n)e2, mei + (L + 2 — n)e2 (region inside the dashed-dotted line in Figure E])). 



Lemma 4.14. Letm,n < L/2. If'qE Vqi^e-nw , F = andVx ' = 

for all X such that {xe2 = L + 1 — n,0 < xei < m) and for all x such that 
[xci = m, L + 1 — n< xe2 < L + 1), then r] E Vr^ . 

Proof. Following the same lines of the proof of Proposition 14.101 it is 
immediate to verify that if and K^^ are empty, then 821- Lei U (e2 + 

qNE-nw-^ is emptied in at most L steps. Since the hypothesis exclude the 
existence of a NE (SE) path for (K^^) of length larger or equal than 

L, the same arguments used in Lemma [4.111 allow to conclude that both this 
sites are emptied in (at most) steps and the proof is concluded. Note that 
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Figure 7: The regions TZl and its partition in the grid of (? squares, with 
c = L/l (here c = 6). The region T^.^ is depicted inside the dashed hne. 
The horizontal region plus the two ending triangles which is delimited by the 
dashed dotted line is Ri (see condition (v) in the text). Instead, the two 
vertical regions delimited by the dashed dotted line form Ci (see condition 
(iv) and (vi) in the text). 

all the hypothesis do not involve the value of the occupation variables inside 

Remark 4.15. Analogous sufficient conditions in order to expand the bot- 
tom (right or left) half of Ol, towards the bottom (right or left direction, 
respectively) can be established by applying the invariance of constraints un- 
der rotations of 90 degrees. 

We are now ready to prove Lemma 14.71 

Proof of Lemma \4.7\ Choose two integers i and L such that L = ci with 
c > 1 also integer. Then divide TZl into a grid of c columns and c rows 
of width i. Let Qi, ■ ■ ■ Qc'^ be the squares identified by this grid and inside 
each square consider a centered octagon Or with r < i/2 and let be the 
octagon associated to Qi (octagon inside the continuous line in Figure [8]^a). 

We will now state a series of requirements (i)-(viii) which involve only 
the configuration inside TZl and are sufficient in order to guarantee that we 
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Figure 8: (a) Zoom on square Qj: The empty internal octagon Ol (guaranteed 
by condition (i)) and (inside the dashed hne) the empty region O^j^ reached 
thanks to conditions (ii) and (iii).(b) Region inside the bold continuous line 
is emptied thanks to the conditions (i)-(vii) 

can empty TZl in (at most) steps. The sufficiency of these conditions can 
be directly proven by using step by step the previously proved properties (as 
detailed below) 

(i) There exists at lest one Qi with c < i < — c such that O). is 
completely empty. Denote by / be the smallest integer such this property 
holds and let K^^ , K^^, K^^ and K^^ be the NW, NE, SW and SE key 
sites corresponding to the octagon Ol- Note that the conditions on i exclude 
that Qi belongs to the top or bottom row of squares. 

(ii) The configuration in Qi belongs to the events /C^^, /C^^,/C^^ and 
/Cf^. This, as has been proved in Lemma [4. IH is sufficient in order to expand 

of one step, i.e. to reach an empty octagon Or+i! 

(iii) The configuration in Qj verifies also the constraints (Lemma 14.111) 
required to expand further the empty octagon until reaching the border of 
Qi. Let 0^12 be the empty region reached via this procedure (region inside 
the dashed line in figure [8] (b). 

(iv) Consider the left and right half of Olj^- The sufficient conditions 
of Lemma 14.141 (rotated of 90 and 270 degrees, see Remark 14.151) with m = 
n = r hold in order to expand these regions (towards the left and the right, 
respectively) via subsequent one site steps until emptying all sites in the 
region Rj composed by a row of squares plus two triangular region at the 
right and left which are at distance i/2 from the border of TZl parallel to 
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ei + 62 (-R/ is the horizontal region dehmited by the dashed-dotted hne in 
Figure [7]). Note that at each step this sufficient conditions does not involve 
any of the occupation variables inside for j ^ I. 

(v) Let Qj C Ri {Qm C. Ri) be the leftmost (rightmost) square such that 
a vertical line drawn from J (M) crosses the bottom (top) border of TZl-, as 
shown in Figure [7] (the existence of such Q j and Q m is guaranteed, since / 
does not belong neither to the top nor to the bottom row). The sufficient 
conditions of Proposition 14.141 and Remark 14.151 with m = n = r holds in 
order to expand with subsequent one site steps the bottom half (top half) of 
^//2 until touching the bottom (top) border of TZl- Note that again 
at each step these conditions do not involve any of the occupation variables 
inside for j 7^ /. Let Cj be the region emptied via this procedure (which 
is composed by the two vertical regions delimited by the dashed-dotted line 
in Figure [7]). 

(vi) For each square Qj such that Qj fl C/ 7^ we repeat the same 
procedure applied in (iv) to empty an horizontal region with the same shape 
of Rj by imposing necessary conditions which do not involve any of the 
occupation variables inside for j 7^ /. 

In Figure [U^b) we depict (inside the bold continuous line) the region which 
is emptied thanks to all the previous requirements. This includes the region 
TZl \ {Ri U Rr) (region inside the dashed line of Figure [8]), where Ri and Rr 
are two lateral strips of width £3/2 (regions delimited by vertexes A,E,H,D 
and F,B,C,G in Figure [H]). 

Let 

Ai:={x:x = -(L-3/2£+z)ei-L/2e2+a(ei+e2), with (i-l)Lc/6j <a< i[c/6\) 

Bi:={x:x = -(L-3/2£+z)ei-L/2e2 + a(ei + e2), with z[c/6j +1 < a < L) 
Di ■={x:x = (L-3/2£+i)ei+L/2e2-a(ei+e2),with (i-l)[c/6j < a < i[c/Q\) 
Ei:={x:x= (L-3/2^ + z)ei +L/2e2 -a(ei + 62), with z[c/6j + 1 < a < L) 

Note that ujf/(v4i U B,) = Ri and uJf/(A U Ei) = Rr. 

(vii) There exists at least one empty site inside each Ai for all i G [1,3/2^]. 
This implies that we can for sure empty all sites in U'^^^Bi. The proof can 
be immediately obtained by subsequent applications of Proposition 14.51 

(viii) There exists at least one empty site inside each for all i G 
[l,3/2£]. This implies that we can for sure empty all sites in U^/^^Ei. The 
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proof can be immediately obtained by subsequent applications of Remark 

The proof of the sufficiency of conditions (i)-(viii) in order to empty TZl is 
then completed by noticing that the union between the region TZl \ {Ri U Rr) 
(emptied via conditions (i)-(vi)) and the regions U^^^Bi and U^^^Ei (emptied 
via conditions (vii) and (viii), respectively) covers TZ^, namely TZ^ C (T^l \ 
{RiURr))ut\' {B^UEi). 

In order to complete the proof, we are now left with evaluating the prob- 
ability of such conditions. If we now denote by P{j) the probability (w.r.t. 
/i^) that the property stated in point (j) is satisfied we get 

e-r 

P{ii n in) > JJ 
i=i 

P{w)P{v)P{vi) > ] 

P (vii) = P (viii) = (1 -pLc/6j)3/2/ 

where the second and third bound follow using Lemma [4.131 and Lemma [4. 141 
and for simplicity of notation here and in the following we drop the index 
OP from the oriented percolation correlation length. 

Note that we have chosen conditions (i)-(viii) in order that the event 
defined by (i) is independent from all the others, since it is the only condition 
which involves the configuration on the small octagons centered inside the 
squares of the grid, i.e. the O^s. On the other hand the events required by 
conditions (ii-viii) are positively correlated (they are all non increasing event 
under the partial order rj -< t]' ii ri{x) < rj'{x) Vx). Thus by using again FKG 
inequality and the above inequalities we get 

l^f{{Tlj^'r] e V^J = I-iP{Gl) > P(i)P{iifMii)P{iv)P(y)P(yi)P{yiif 

(4.39) 

Let r := 2{log^ > ^, c := (1 - p)-35'i°g«', £ := ^Mog^ and define the 
function L(p) as follows 

L(p):=(l-p)-3«^'°s«Vloge (4.40) 



exp( 



exp( 
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From the above inequality we get in leading order as ^ — > cxo 

P(z) > 1 -exp(-(l-p)^2(«i°s«)') 

P{tt) > exp(-4/0 
P{w)P{v)P{vi) > 1 - exp(-^3logO 
P{vii) = 1 -eMogeexp(-|logp|/4(l 

It is immediate to verify that in the limit p p^^, since ^ — > cx) all these 
quantities go to one and for any e > there exists p(e) (with p(e) 
when e \ 0) such that ^ := ^(p(e)) is sufficiently large to guarantee that the 
product on the left hand side of fl4.39p is bounded from below by 1 — e for 
P > p(e) • This implies the result of Lemma 14.71 with 

L{p,e)=L{p) (4.41) 

for p > p{e), where L{p) has been defined in fl4.4ip . The result for all the 
densities, p < p(e) < p'^^, and with the choice L{p,e) = L{p{e)) trivially 
follows from attractiveness which implies monotonicity of the probability of 
the good event (on a fixed size) under p. 

5 Discontinuity of transition: proof of Theo- 
rem 13.41 

The proof is composed of two steps. First we construct a set of configurations 
such that the origin is frozen, i.e. it cannot be emptied at any finite time. 
Then we prove that this set has finite probability at p = p^^- This is the same 
strategy which we already used to prove the upper bound of Pc- However the 
clusters which block the origin will now be of different type, the key feature 
being the existence of North-East and North- West occupied clusters which 
are mutually blocked. It is thanks to these structures that the properties of 
the transition are completely different from those of oriented percolation. On 
the other hand, we will see that an important ingredient which guarantees 
a finite weight to these configurations is the anisotropy of typical blocked 
clusters in each one of the two directions, i.e. anisotropy of conventional 
oriented percolation (it is here that Conjecture 13.21 is used). 

Before entering the proof of Theorem 13.41 let us establish a result which 
will be used here as well as in the proof of Theorem 13.51 (ii). Recall that we 
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denote by ^ a rectangular region with sides of length a parallel to ei + 62 
and sides of length h parallel to — ei +62- Let Ri be a rectangular region A^ ;, 
centered at the origin and let X := xi, . . . , a;„ with Xi G -Ri be a North-East 
path which spans -Ri, i.e. with xi and Xn belonging to the two sides of -Ri 
which are parallel to — ei + 62 (see FigjO]). Choose c and d such that 2c < a 
and d> h and let R2 and -R3 be two rectangular regions of the form A^ which 
are centered on the line ei + 62 with centers at (a — c )/(2y2)(ei + 62) and 
— (a — c)/(2-\/2)(ei + 62) (we consider without loss of generality that a, h, c, d 
are integer multiples of 2a/2). Note that both R2 and R^ intersect Ri and 
both sides of Ri which are parallel to — ei + 62 lie on the two more far apart 
sides of i?2 and R^ along this direction. Finally, let Y := {yi, . . . ,1/^) and 
Z := {zi, . . . , Zm') with tji G -R2 and Zi G R3 be two North- West paths which 
span i?2 and -R3 respectively. Lemma [5?T] below states that if the three above 
defined paths X, Y and Z are occupied, the subset of X which belongs to 
Ri \ {R2 U -R3) cannot be erased before erasing at least one site which belongs 
either to Y or to Z (actually along the same lines of the proof below one can 
prove the stronger result that one needs to erase at least M := {d-b)/{2y/2) 
sites in either Y or Z). Let X := X \ {R2 U -R3) (sites inside the dashed 
contour in Figl9]) and ta be the first time at which at least one site in A is 
empty. The following holds 

Lemma 5.1. If rio{w) = 1 for all w e {X UY U Z), Tx > Tyuz + 1 
Proof 

In this proof for simplicity of notation r-^ will be sometimes simply denoted by 
r. Let / and u to be the indexes such that xi G i?i\i?3, xi-i G -R3, Xu G Ri\R2 
and Xu+i G -R2 (it is immediate to verify that {u — I) > (a — 2c)/\/2). The 
definition of implies that r]r-i{xi+i) = ?7r~i(xj+i) = 1 for / < z < n, thus 
Vri^i) = 1. This in turn implies that r]r{xu) = or rir{xi) = (or both). We 
will now prove that if the first possibility occurs this implies that at least 
one site of Y should be emptied before r. The other case, r]r{xi) = 0, can be 
treated analogously leading to the result that at least one site of Z should 
be emptied before r. In both cases the result of the Lemma follows. The 
first observation is that, since we consider the case r]r{xu) = 0, this implies 
rir-i{xu+i) = 0. This follows from the fact that rir-i{xu-i) = 1, therefore 
the SW neighbours of Xi are not both empty at r — 1. Thus both its NE 
(which include Xu+i by definition) should be empty at r — 1, otherwise the 
emptying of x„ at time r could not occur. This procedure can be iterated to 
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Figure 9: The three rectangles -R2 and i?3 with the corresponding occu- 
pied clusters, X, F and Z. The red (blue) line correspond to and iy and 
their intersection point is P. Here P ^ X, P , therefore we are in case 
(b) for the proof. Black dots inside the dashed contour belong to X, i.e. the 
set of sites which are guaranteed to be occupied before at least one site in 
Y VJ Z has been emptied. 

prove that if Su+i is the first time at which Xu+i is emptied it verifies 

Su+i < Su+i-i - l<--<s„-l = r- l (5.42) 

for < 2 < n — M. Since both X and Y span R := RiCl R2 and X connects 
the two sides of R which are parallel to —61+62 and Y those that are parallel 
to 61 + 62, if we denote by ix (^y) the continuous line obtained by joining 
the sites in X (in Y), it is immediate to verify that £x and iy do intersect. 
Denote by P an intersection point for ix and iy. Since ix is by construction 
composed only by segments of the form 62 and 61+62, while iy is composed 
by segments of the form —61 and —61 + 62 it can be easily verified that (a) 
either P belongs to the lattice (and therefore to both X and Y) (b) or it does 
not belong to the lattice but P ± (61 + 62)/2 belong to X and P ± (61 — 62) 
belong to Y (this is for example the case for the paths depicted in FiglH])- 
Let us treat ) and (b) separately. 

(a) Since P G X and P G P2, we can identify an index j such that 
u < j < n and P = Xj. Therefore by using fl5.42p . we get that the first 
time at which P is emptied, s{P), verifies s{P) = sj < Tj^ — 1. Since P also 
belongs to Y we have ryyz < Sj and therefore t-^ > Tyuz + 1- 

(b) Since P — (6i + 62)/2 G X and P G (P2UXu), we can identify an index 
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j such that u < j < n and P — {61 + 62}/ 2 = Xj. By using fl5.42p . we get that 
Sj < — 1 and Sj_i > Sj. Therefore both the NE neighbours of Xj should be 
empty at time sj — 1 (since its SW neighbours are not both empty). These 
NE neighbours include P + (ei — 62) /2, which therefore would be empty at 
time Sj — 1. Since P + (ei — e2)/2 G Y, we have proven that tyuz < — 1. 
Putting above results together we conclude that t-^ > Sj + 1 > Tyuz + 2. 

□ 

We can now proceed to prove the discontinuity of the transition. 

Proof of Theorem 13.41 Fix = ii > and let Aq be the rectangle with 
the shape of A^^ i/i2£o centered at the origin and A\ and A^ be two rectangles 
with the shape of A^j/12/i centered on the line 61 + 62 at distance £0 — l/24£o 
from the origin, as shown in Figure (TUl Let also ii = 2£j_2 and consider 
two infinite sequences of the rectangular regions A^^^ and Af with the shape 
of Ae^^i/uii for i even and of Ai/124,4 for i odd. As shown in Figure [TOl the 
centers of A;^ for i odd (even) lie all on the Ci + 62 (— Ci + 62) line and their 
distance form the origin satisfies |c4| = ii-i/2 — ti/2A. 

Let Si for i even (for i odd) be the event that there exists an occupied 
NE [NW) path which connects the two short sides of A^. Let also Sq be the 
event that there exists an occupied NE path which contains the origin and 
connects the two short sides of Aq. 
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Figure 10: We draw the first elements of the two sequences of increasing 

1 2 

rectangles TZ^ . The figure is rotated of 45 degrees for sake of space, the 
coordinate directions ei and 62 are indicated. 

The origin is frozen if r/ G 5*0 S} S^, namely 

Poo{p)>f^'{Sor}Z,slr}T=iSl) (5.43) 

Indeed if we let Al be the subset of the path spanning A^ which belongs 
to A^ \ {A\j^^ U A|_,_j^) and Xj be the first time at which at least one of the 
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sites in A} or is emptied, by using Lemma 15.11 (and the analogous result 
for the structure rotated of 90 degrees) it can be easily established that: (i) 
Ti > Tj+i + 1; (ii) the origin cannot be emptied before time ri. For any 
r > if we choose i > t, from (i) (and the fact that Tj > 1) we conclude 
immediately that ti > i > t. Therefore from (ii) the result T^?7(0) = 1 
immediately follows. This, together with the arbitrariness of r immediately 
leads to fl5.43p . Therefore, by using FKG inequality and definition we get 

Poo(p) > ^^(5o) n (5-44) 

i=l,...oo 

By using the same mapping of NE oriented paths to paths of oriented perco- 
lation used for Lemma [4.11 (and an analogous version in the NW direction), 
it is then immediate to verify that n^{Sl) = fi^{Si), where Si are the oriented 
percolation events defined in fl3.18p . Therefore the result Poo{Pc) > follows 
from Conjecture 13.21 □ 
Note that the structure which we have used to block the origin does 
not contain any infinite cluster neither in the North-East nor in the South- 
West directions, therefore it could be emptied if we had chosen to block 
just in one direction (e.g. by taking Ax '■= {t] : r] E {S^^ U S^^)) or 
Ax := {rj : T] E {S^^ U £^^))). Oriented percolation corresponds to this 
choice of blocking just in one direction, this is why our proof of discontinuity 
does not apply to this case where it is indeed well known that the transition 
is continuous. 



6 Finite size effects: proof of Theorem 13.5 



In this Section we will prove Theorem 13.51 which provides upper and lower 
bounds for the scaling of the crossover length H(|o) (see definition fl2.12p ) in 
terms of the correlation length of oriented percolation. 
Proof of Theorem 13.51 (i) 

Let A2L and Al/2 be two squares centered around the origin and of linear 
size 2L and L/2, respectively. If we recall the definitions given in section 
14.21 for regions Rl and Rl, it is immediate to verify that Al/2 C Rl and 
Rl C A2L. This implies that if we take a configuration 77 in A2L and we 
evolve it with occupied boundary conditions, Al/2 is completely empty if the 
restriction of rj to Rl is a good configuration. Thus E{L,p) > pP{G^) and 
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the upper bound for the crossover length follows immediately from Lemma 
Oand equation (531]) • 

Proof of Theorem 13.51 (ii) 

Let s and n be two positive integers and consider a square centered around 
the origin of linear size 4ns with two sides parallel to ei + 62 and the others 
parallel to — ei + 62, namely a region Ains,ins- Inside this square we draw 
a renormalized lattice with sides parallel to the square sides and minimal 
step 4s as shown in Fig. [11] and, without lack of generality, we let the origin 
belong to this renormalized lattice. Then we draw around the origin the 
structure shown in Fig. [11] (a) which is composed by the intersection of four 
rectangles, two of the form Ag^ (dashed-dotted line in Fig. [TTT a) and two of 
the form (dashed line in Fig. [TT] (a). Analogously, around the four sites 
n.n to the origin, we draw the structure obtained by a reflection from the one 
around the origin (so that for each of these sites one of the four rectangles 
coincides with the one of the origin). Then we continue this procedure until 
depicting one such structure centered on each site of the renormalized lattice 
in such a way that the structures on two neighbouring sites coincide always 
up to a reflection. This leads to the final structure depicted in Fig. [TT] (b) 
which contains 2N = 0(2n^) rectangles: half are of the form Aqs,s and will 
be denoted by 7?.i, . . . TZn with TZi being one of rectangles around the origin; 
and half are of the form A^ g^ and will be denoted by T^-at+i, . . .TZ2N- For 
1 < i < N {i > N), we let Si be the event that there exists an occupied 
NE (NW) path which connects the two short sides of TZi. If we evolve the 
dynamics with occupied boundary conditions on the square A2^„^ inscribed 
inside A^n.s,4ns (square inside the dotted line in Fig. [TT](b)). the inner square 
^V2/2ns ^2V2ns '^^^ never be completely emptied (provided n is sufficiently 
large, n > 12, in order that A^^^ns contains at least one rectangle Ri). In 
other words, recalling definition fl2.11l) . the following holds 

E{V2ns, p) < 1 - fi'{n^,=iSi) (6.45) 

The key observation to prove fl6.45p is the following. Choose a rectangle TZi 
and focus on the four perpendicular rectangles by which it is intersected. Let 
71 j and TZ'^ be the two that are more far apart (see Fig. [TT] (b)). li i < N 
{i > N), let Ai be the subset of the NE (NW) spanning cluster for TZi which 
is contained in TZi \ ("^j U7^^ ). The union over all rectangles of the sets Ai has 
the property that none of these sites can ever be emptied. This is an easy 
consequence of the result in Lemma 15.11 (and of having imposed occupied 
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Figure 11: (a) A^ns,4ns (here n = 4) and the structure around the ori- 
gin. (b)The dashed and dotted-dashed rectangles form the structure described 
in the text obtained by drawing around each renormahzed site either the same 
structure as for the origin or the one obtained by a reflection symmetry. The 
dotted square corresponds to f^2V2ns- 

boundary conditions on Ag^^^). Therefore since A^^^ns ^ {^1=1^1) ^ 
(provided n > 12), fl6.45p immediately follows. 

Again we can use FKG inequality and conclude that 

E{V2ns, p) < 1 - fi^iSiY^ (6.46) 

Note that now all rectangles have the same size, namely fJ^^^Si) does not 
depend on i. The probability that there does not exist a NE cluster spanning 
TZi is clearly bounded from above by the product of the probabilities that 
on each of the s^~^ slices of the form Aqs,s^ which compose TZi there does 
not exist such a cluster when empty boundary conditions are imposed on 
the long sides of the slice. Therefore if we recall Conjecture 13.11 and we let 
■s = ^(p), we get 

lim 1 - fiP{Si) < lim (1 - c^opf'^ 

P/Pc P/Pa 

By using this inequality inside fl6.46p we immediately get 

lim E(V2<(p),p) < lim 1 - exp \-^{n + 1)^ exp(-a^^~^)l 

P/Pc P/Pc 
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(6.47) 
(6.48) 



where a = \ log(l — Cqp)\ (and we used the fact that (1 — c^p)^^ " ^ and 
log(l - x) > -2x for < X < 1/2 and 2N < 2{n + 1)^). Therefore, if we let 
n{p) = l/Ae:icg{a^{pY~^ /2) we get 

hm E(ci^exp(c2^^"^),p) < 1/2 (6.49) 

P/Pc 

with ci = 1/(2a/2) and C2 = | log(l — c^p) |/2. By recalhng the definition of 
the crossover length fl2.12p the proof is concluded. □ 
As discussed in the Section 2 another possible definition of the crossover 
length would have corresponded to defining E{L,p) in 02.121) as the proba- 
bility that the origin is empty in the stationary configuration which is reached 
after steps when we evolve the configuration with filled boundary condi- 
tions on Al. In this case the previous proof would have to be modified. First 
one has to traslate the structure in order that the origin is at the center of 
rectangle Ri. Then the events Si for i > 2 remain unchanged but Si should 
now be the event that there exist a structure which contains the origin similar 
to the one used in the previous Section to prove discontinuity (see Fig. [TD]) 
and this up to the size of Ri (i.e. with the longest rectangles Aj spanning Ri 
in the parallel direction). Then by combining the arguments used to prove 
Theorem 13.41 and the above Therem 13.51 (ii) it is not difficult to see that the 
origin is blocked by this structure, thus leading again to inequality fl6.45p . 
Finally, to establish the upper bound fl6.49p for the new E{L,p), one needs 
to use a conjecture on the properties of oriented percolation which is slightly 
more general than 3.2 and is again expected to be correct on the basis of 
finite size scaling and numerical simulations. 

7 A related Kinetically Constrained Spin Model 

In this section we define a Kinetically Constrained Spin Model that we 
introduced in [17], which is related to the cellular automaton that we have 
considered in previous sections. Configurations 77 are again sets of occupation 
variables rj^ G {0, 1} for x ^1? distributed at time with /x'', but evolution 
is not deterministic. Dynamics is given by a continuous time Markov process 
with generator L acting on local functions / : ^ M as 

I^J{ri) = Y,cM\f{rr)-f{ri)\ (7.50) 
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with 

V:-=i ^'"^^ '^^Z^ (7-51) 
I r]z n z X ] 

The rates c^(?7) are such that the flip in x can occur only if the configuration 
satisfies the same constraint that we required for the cellular automaton in 
order to empty the same site, namely 

( if r/^ ^ A 
Cxiv) ■= { P if ?7 e A and r]^ = (7.52) 

1 — p a T] G Ax and rjx = I 

It is immediate to check that the process satisfies detailed balance with re- 
spect to fi'^, which is therefore a stationary measure for the process. This 
property is the same as for the process without constraints, namely the case 
Cxiv) = p(l ~ vi^)) + (1 ~ P)l{^)i but important differences occur due to 
the presence of constraints. In particular for our model /x^ is not the unique 
invariant measure. For example, since there exist configurations which are 
invariant under dynamics, any measure concentrated on such configurations 
is invariant too. A direct relation can be immediately established with the 
cellular automaton studied in previous sections: configurations which are left 
invariant by the stochastic evolution are all the possible final configurations 
under the deterministic cellular automaton evolution (since all sites in such 
configurations are either empty or such that if ri{x) = 1 than rj ^ Ax)- A 
natural issue is whether on the infinite lattice, despite the existence of sev- 
eral invariant measures and of blocked configurations, the long time limit of 
all correlation functions under the Markov process approaches those of the 
Bernoulli product measure for almost all initial conditions. By the spectral 
theorem this occurs if and only if zero is a simple eigenvalue of the generator 
of the dynamics, i.e. if Cfo{r]) = with /o G L2{np) implies that /o is con- 
stant on almost all configurations, i.e. on all except possibly a set of measure 
zero w.r.t. /i'' (see [12] and, for a specific discussion on Markov processes 
with kinetic constraints, see Theorem 2.3 and Proposition 2.5 of [3]). 
By the result in Lemma 14.11 it is immediate to conclude that the process is 
for sure not ergodic for p > p'^^, since for example the characteristic func- 
tion of set j^^^-^'^ ig an invariant function which is not constant. On the 
other hand. Lemma 14.21 establishes irreducibility for p < p^^, namely that 
a.s. in pp there exists Va; a path rji = r], . . .rjn = t]^ such that r/j+i = rjf and 
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Cy{f]i) = 1 B Ergodicity for p < p^^ can then be immediately established 
thanks to irreducibility and the product form of Bernoulli measure (see [3] 
Proposition 2.5). Therefore at p^^ an ergodicity breaking transition occurs. 
Furthermore, by using the result of our Lemma W7j[ in [4] it has been proved 
that the spectral gap of this process is strictly positive at any p < p^^, i.e. 
correlations decay exponentially in time to their equilibrium value. 
Finally, for the connection with the physics of liquid/glass transition, let us 
analyze the Edwards- Anderson order parameter q which corresponds to the 
long time limit of the connected spin-spin correlation function. 



g = V fi^'ivhiv) = Y] f^^iv) lim lim 

— ^ — ^ t—*oo |A|^oo 

V V 



J2\M Mx)r]t{x)-p^ 



(7.53) 

where denotes the mean over the Markov process started at rjo = r). 
Ergodicity guarantees that q = for p < p^^ . For p > p^^ one can obtain 
for a fixed initial configuration r] the following inequality 

q{ri) > {1- p)'^pob/il- pob) 

where p is the fraction of occupied sites in t] and pob is the fraction of oc- 
cupied sites remained after performing on rj the emptying process defined 
by the cellular automaton. Under the hypothesis that both these quantities 
have vanishing fluctuations in the large L limit with respect to the Bernoulli 
measure for the initial configuration, one finds that for p > it holds 

and therefore q{pc) > 0. 

Finally, as explained in [1 7] , the fact that the crossover length for the cel- 
lular automaton diverges faster than exponentially toward the critical den- 
sity should correspond to (at least) an analogous divergence in the relaxation 



■^This follows immediately from the proof of Lemma 14.21 we have shown that for the 
cellular automaton if we consider a sufficiently large finite lattice A l around x there exists 
a path from i] which subsequently empties all sites in Al. For all the moves rji rji+i ~ rjf 
in this path the rate cf is non zero since the constraint for the cellular automaton and the 
stochastic process are the same. We can construct analogously a path which goes form rf 
to a configuration which is empty in A^. Then, by connecting the two paths we get an 
allowed (i.e. with strictly positive rates) path from rj to rj^ . 
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times for this stochastic modeL 

As discussed in [17], the first order/critical character of this dynamical tran- 
sition is similar to the character experimentally detected for liquids/glass 
and more general jamming transitions. To our knowledge, this is the first 
example of a finite dimensional system with no quenched disorder with such 
a dynamical transition. 
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